Optical angular momentum in dispersive media 
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The angular momentum density and flux of light in a dispersive, rotationally symmetric medium are derived 
from Noether's theorem. Optical angular momentum in a dispersive medium has no simple relation to optical 
linear momentum, even if the medium is homogeneous. A circularly polarized monochromatic beam in a 
homogeneous, dispersive medium carries a spin angular momentum of ±h per energy huj, as in vacuum. This 
result demonstrates the non-trivial interplay of dispersive contributions to optical angular momentum and energy. 
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The discovery by Allen et al. [1] that light beams can carry 
orbital angular momentum has led to intensive study of the 
angular momentum of light and its applications (see yfl for 
a recent review). In vacuum, the linear momentum density 
p = £oE x B of light determines its angular momentum den- 
sity r x p. A beam with azimuthal phase dependence car- 
ries orbital angular momentum in the direction of propaga- 
tion, distinct from the more familiar spin angular momentum 
of circular polarization (±7j per photon) [2]. Inside a medium, 
however, these expressions for the linear and angular momen- 
tum densities of light are no longer correct. If the beam has 
frequency components for which absorption is significant then 
there is no conserved optical linear or angular momentum. 
But if the beam is confined to a frequency range where losses 
are negligible, then light in a rotationally symmetric medium 
will carry a conserved angular momentum. A frequency range 
with negligible losses will in general exhibit dispersion, visi- 
ble light in glass being the most famous example. The prob- 
lem solved in this paper is the following: what is the optical 
angular momentum density and flux in a dispersive, rotation- 
ally symmetric medium for beams with negligible absorption? 

Even if the beam is monochromatic, dispersion contributes 
to the energy-momentum and angular momentum of light. 
This is best known from the Brillouin expression for the 
time-averaged energy density of monochromatic light in a 
dispersive medium y|], which depends on the derivatives 
de(uj) I dui and dfj,(ui) /dui of the permittivity and permeability 
at the frequency of the beam. The time-averaged monochro- 
matic momentum density also depends on de(ui)/duj and 
d^(uj)/duj flU. 

The energy density and momentum den- 
sity for beams with a finite frequency range has a much more 
complicated dependence on the dispersion ]4J]. Here we find 
the angular momentum density and flux for a finite frequency 
range and also specialize the result to the monochromatic 
case. Unlike in vacuum, the angular momentum density has 
no simple relation to the linear momentum density (even in 
a homogeneous medium where the latter obeys a conserva- 
tion law). This means that the conserved angular momentum 
cannot be constructed from the energy-momentum tensor, al- 
though this is how optical angular momentum in media is usu- 



ally addressed Jj§]. The conserved optical angular momentum 
is associated with rotational symmetry and must therefore be 
calculated using Noether's theorem. 

Note that the problem solved here is self-contained, unam- 
biguous, and unaffected by considerations of angular momen- 
tum transfer between light and matter. We consider only the 
correct expression for conserved optical angular momentum; 
an analysis of the transfer of optical angular momentum to 
matter requires ingredients not needed here BUt]. The ques- 
tion we address is one of the most basic that can be posed for 
any field theory. The solution to the problem, in combination 
with the result for the energy-momentum tensor, demonstrates 
a remarkable interplay of dispersive contributions to the con- 
served quantities of light. 

As in [4], our results are derived from Noether's theorem 
applied to the electromagnetic action in a dispersive medium, 
assuming a finite frequency range with negligible losses. In 
such a frequency range the dielectric functions can be fitted to 
an even series in frequency: 
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£ (r, W ) = ^ £2n (r) W 2 ", K (r, W ) = ]T K2n (r) W 2 ", (1) 

n=0 n=0 

where e(r, w) is the relative permittivity and the relative per- 
meability is jtx(r, lu) — n(r, In practice the series in ([TJ 
will have a finite number of terms and will represent a fit to 
dispersion data in the frequency range of interest, but we al- 
low an infinite number of terms to obtain the results in their 
greatest generality. The expansions (JTJ are standard in treat- 
ing dispersion in frequency ranges where absorption is neg- 
ligible (see e.g. B9D). To ensure the existence of a conserved 
optical angular momentum, the medium is taken to be rota- 
tionally symmetric, with dielectric functions ([TJ depending on 
r = y/r ■ r. The D and H fields in the frequency domain and 
in the time domain are 

D(r, uj) = e Q e(r, w)E(r, u), D(r, t) = e e(r, id t )E(r, t), 
H(r,w) = K n(r,u)B(r,uj), H(r,t) = K n(r, id t )B(r, t). 

where kq = fa . The electromagnetic action in the medium 
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with the dynamical variables taken to be the scalar potential . 
and vector potential A, defined by 



E = -VcA-&A, B = VxA. 



(3) 



Variation of <p and Am© gives the macroscopic Maxwell 
equations (with no free charges or currents) 

e V.[e(r,*9t)E]=Q, (4) 
k V x [re(r, i<9 t )B] = e e(r, id t )d t E. (5) 

The other two Maxwell equations are identities due to (|3). 

The action (O is invariant under active rotations of the dy- 
namical fields and A around the origin r = 0, and this in- 
variance implies the existence of a conserved quantity, optical 
angular momentum. For a homogeneous medium the invari- 
ance of the action holds for rotations of the dynamical fields 
around an arbitrary point, but this point can always be chosen 
as the coordinate origin. An infinitesimal rotation of the scalar 
field (f> around the origin is given by 



4>(r,t) -> <t>(r-6r,t), 



(6) 
(7) 



where fly = —Clji denote the three independent infinitesimal 
parameters of the rotation. (We use tensor notation through- 
out, with indices lowed and raised by the metric and its 
inverse g 13 .) The vector potential undergoes the rotation 



and the infinitesimal variations of <b and A are thus 
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8(j> = -fi^VV, SAi = <>,,,1' - n jk r k V j Ai. (9) 

Noether's theorem guarantees that if we let the rotation 
parameters depend on space and time, i.e. Qjj = f2y(r,i), 
then the change in the action (O under (0 can be written in 
the form 



6S = -- j (fx (L ij d t n i:j + M h, 'YiM 



(10) 



where the angular momentum density L lJ = — U % and flux 
M kij — —M kjl obey the conservation law 



d t L ij + V fe M^ = 0. 



(11) 



The angular momentum density L tJ emerges as an antisym- 
metric tensor; its dual is a vector and will be constructed 
later. Integrations by parts, in which surface terms are to be 
dropped, are required to achieve the form (TTOb . Use must also 
be made of the identities (10) and (1 1) in ref. |4J]. The result is 



as follows, where we use square brackets to denote antisym- 
metrization of indices [ 14], is the completely antisymmet- 
ric Levi-Civita tensor, and C is the Lagrangian density (i.e. the 
integrand in the action (0): 



L i:j =2e A [j e(r,id t )E i] + 2e Qr[i W j] A k e(r, id t )E k 
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The angular momentum density ( TT2l and flux (13[ satisfy the 
conservation law (fTTT) when Maxwell's Eqs. (01 and <(5j are 
used; one obtains 



d t Ui + V fe M^' =e Q r^E k \v j h(r,id t ) 
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which vanishes because A l V^f(r) — for any function of 
r = y/r ■ r. 

The results (TT2l and (fT3l are not gauge invariant. Equiva- 
lent gauge-invariant expressions are found by noting that the 
quantities 

F kij = 2e r^A^e(r,id t )E k , (15) 
G lkl i = AK r lj A^n{r, id t )V [k A^ (16) 

identically satisfy 

d t V k F kl3 + V k {-d t F klJ + ViG lklJ ) = 0. (17) 

From (TTTb we see that adding V k F kl; > to L tJ , and adding 
_Q t pkij + Yjgrifeij to M ktJ , does not affect the conservation 
law (fTTT) . Performing these additions, and using Maxwell's 
Eqs. © and ©, we obtain the following gauge-invariant an- 
gular momentum density and flux: 



Lij =2r [iPj] 



5 ° E E (-l) n+m ean W- 1 ^"-" 1 ^ 



n—l rn — 1 

oo 2n 



«oEE (-l) n+TO «a»(r)a? l - 1 B [< ^ n - ,n B i ] > 

(18) 



n—l rn—1 



3 



M; 



klj 



(19) 



where pi and aij are, respectively, the linear momentum den- 
sity and stress tensor in a dispersive medium 0]: 

Pi = (D x B)i 
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(20) 



ofc = - EiDj - HtBj + - gij (E D + B H) . (21) 

The linear momentum density d20b is not the Minkowskii mo- 
mentum D x B, which is only valid if there is no dispersion. 
The stress tensor (ETb has the same expression as in the non- 
dispersive case, but is not symmetric in the time domain. The 
momentum density d20l > and stress tensor ( f2TT > obey a conser- 
vation law only if the medium is homogeneous [4]. A striking 
property of the angular momentum density (TT~8T > is that it has 
no simple relation to the linear momentum density p, even 
when linear momentum is conserved. The angular momen- 
tum flux ([T9l >. on the other hand, is completely determined by 
the linear momentum flux, the stress tensor uij. The result 
(|T~9T > is the same expression for the angular momentum flux 
as in vacuum Q], except that the stress tensor is symmetric in 
vacuum so that rpayjj. = cr^ra, whereas this last equality 
does not hold in a dispersive medium. 

It is more usual to write optical angular momentum den- 
sity as a vector, with the flux as a second-rank tensor; these 
quantities are simply the duals of and M k ij , defined by 



V = le l3k L 3k , M ij = \t lkl M\ 
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The angular momentum density vector L and second-rank flux 
M' lJ are, from (18}, © and (32}, 
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and they satisfy the conservation law 

<),i: + WjM ji = o. 



(23) 
(24) 



(25) 



The results (TS} and (19}, or their duals (|23]l and d24ji, are the 
solution to the problem posed at the beginning of this paper. 



Similar to the energy density and momentum density in a 
dispersive medium Jj], the angular momentum density sim- 
plifies considerably for time-averaged monochromatic waves. 
We insert the monochromatic E field 
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E(r,t) = -(E (r)e 



c.c 



(26) 



and a B field of the same form into ( 1231 . After time averaging, 
all i-dependent terms vanish and we obtain the time-averaged 
monochromatic angular momentum density L mono : 
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where p mono is the time-averaged monochromatic linear mo- 
mentum density Jj, 0] 



Pmono =Y£(r,o;)Re(Eo x Bq) 



e de(r,u) 



Kq dK(r,u) 



(28) 



The time-averaged monochromatic energy density p mono in a 
dispersive medium is given by the well-known Brillouin for- 
mula nn 
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H ■ Hq. 



(29) 

This energy result is exact only for monochromatic waves, but 
is often used as an approximation in the quasi-monochromatic 
case. The exact energy density for finite frequency ranges 
where absorption is negligible is given in Jj]. Note that the 
energy density expression ( f29b is in fact valid for any spa- 
tial dependence of the dielectric functions, but the momen- 
tum density d28l i only measures a conserved quantity if the 
medium is homogeneous. 

To study a specific example, we consider a homogeneous 
medium, where the dielectric functions e(ui) and k(uj) are in- 
dependent of position. The simplest monochromatic wave in 
this case is a plane wave, but a well-known oddity of a circu- 
larly polarized plane wave in vacuum is that its angular mo- 
mentum density eor x (E x B) in the propagation direction 
turns out to be zero J2l- If the circularly polarized wave is 
given a transverse intensity that falls off to zero at infinity, the 
angular momentum in the propagation direction is found to be 
±h per photon, the expected spin angular momentum yj. To 
capture the full spin angular momentum of a circularly polar- 
ized monochromatic wave we must therefore consider a non- 
trivial transverse profile. The simplest case is to assume that 
the transverse profile decreases to zero slowly compared to the 
wavelength scale yfl. The resulting spatial dependence Eo(r) 
for the monochromatic electric field d26l i. and corresponding 
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B (r), are then [3] 
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We have taken the wave to propagate in the x-direction, with 
{e^, e y ,e z } unit vectors along the Cartesian axes. The trans- 
verse beam profile is given by £(y, z), n p (uj) = y/e(cd)[i(u)) 
is the phase index, k = n p (ui)oj/c is the wave vector and 
the upper/lower signs are for left/right circular polarization. 
One can verify that ( T3Qb and (l3TT l satisfy Maxwell's equa- 
tions when k£(y, z) is much larger than derivatives of £(y, z). 
With the same assumption of a slowly varying £(y, z), the x- 
component of the time-averaged angular momentum density 
d27l for the beam ^30j> and (fJB is 



L x = T£o£ ( w )^U ( rx • V ± )f ± e 01 



/ e{u) dn p (lu) ^ 2 
(32) 

where = {y, z} is the transverse position vector and Vx = 
{V y , V^} is the transverse gradient operator. By integrating 
(l32l over the yz-plane we find the spin angular momentum 
per unit length of the beam. The first term in d32l can be re- 
written using 2£(r^ ■ V ±)£ = Vj_ ■ (rj_£ 2 ) — £ 2 V ±_ ■ and 
■ rx = 2, so the integration gives 



d r±L x = ± 



e e(uj)ng(uj) 
Lun p (uj) 



d 2 r ± £ 2 , 



(33) 



where n g (u>) — d[ujn p (u!))/du> is the group index. The time- 
averaged energy density d29b of the beam (f30b and OTb is 



(34) 



Integration of ( f34t over the yz-plane gives the energy per unit 
length, and dividing this into d33l we find a spin angular mo- 
mentum per unit energy of ±w _1 , just as in vacuum [3J. An 
energy of hcu in the beam thus carries a spin angular momen- 
tum of ±h. 

The time-averaged momentum density (|28l l of the beam 
( f30b and (|3TT i in the x-direction is 
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and dividing this by the time-averaged energy density d34l 
gives a momentum per unit energy of n p {uS)jc. An energy 
of ftuj in the beam therefore carries a linear momentum of hk 
(this is also the result for a plane wave). 

Although we have calculated spin angular momentum and 
linear momentum for an energy Huj, the results presented here 
are purely classical. Nevertheless, we have shown a highly 



non-trivial interplay of dispersive contributions to the classi- 
cal energy, momentum and angular momentum of light that 
indicates the expected results for photon momentum and an- 
gular momentum. A rigorous quantization of light in disper- 
sive media requires a full account of absorption consistent 
with Kramer- Kronig relations [llj]. The energy-momentum 
tensor of classical and quantum light in an arbitrary medium 
obeying Kramer-Kronig relations is given in 111 211 in a form 
that includes the energy-momentum absorbed by the medium. 
The total angular momentum density and flux in an medium 
obeying Kramer-Kronig relations can easily be found from the 
action in 111 UN 211 using Noether's theorem. Separation of total 
energy-momentum and angular momentum into electromag- 
netic and absorbed parts is a difficult problem in general, but 
for frequencies with negligible absorption the results derived 
here and in [4] must follow from a more general treatment that 
includes losses. 

As we stressed in the opening paragraphs, we have not an- 
alyzed angular momentum transfer from light to matter. This 
transfer can occur either by absorption of a pulse carrying an- 
gular momentum or by transmission of the pulse through the 
medium. The results derived here can be used to calculate the 
total angular momentum carried by the pulse when it is com- 
pletely contained inside the medium. An analysis of the angu- 
lar momentum transfered to the material [7] leads to the same 
subtleties that have been long debated in the case of linear 
momentum (the Abraham-Minkowskii controversy fl l"l3"ll ). 
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